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Abstract 

We propose that, within the standard model, the correlation between the ti forward-backward asym- 
metry A t i and the corresponding lepton-based asymmetry Ai - at the differential level - is strong 
and rather clean both theoretically and experimentally. Hence a combined measurement of the two 
distributions as a function of the lepton pt, a direct and experimentally clean observable, would lead 
to a potentially unbiased and normalization-free test of the standard model prediction. To check the 
robustness of our proposal we study how the correlation is affected by mis-measurement of the ti 
system transverse momenta, acceptance cuts, scale dependence and compare the results of MCfm, 
Powheg (with & without Pythia showering), and Sherpa's CSshower in first-emission mode. We 
find that the shape of the relative differential distribution Ai (p l T ) [A t j (p l T )] is only moderately dis- 
torted hence supporting the usefulness of our proposal. Beyond the first emission, we find that the 
correlation is not accurately captured by lowest-order treatment. We also briefly consider other dif- 
ferential variables such as the system transverse mass and the canonical ti invariant mass. Finally, 
we study new physics scenarios where the correlation is significantly distorted and therefore can be 
more readily constrained or discovered using our method. 
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1 Introduction 



Within the standard model (SM), the ti forward-backward asymmetry, A t j, is an interesting variable 
because it tells us about QCD interactions beyond leading order but in a region that should be well de- 
scribed by perturbation theory [1,2]. Furthermore, as the standard model contributions are expected to be 
small [1-5], the measurement of A t i is sensitive to beyond-the-SM (BSM) contributions. The asymmetry 
is quite an unique observable since shifting it requires new physics with non-standard couplings both to 
the tt quark current as well as to the to the current of uu (or possibly dd) initial-state quarks. 1 

The current status of top quark asymmetry related measurements at the Tevatron is intriguing. It is use- 
ful to classify the current data into measurements that directly probe the ti asymmetries and measurements 
that probe daughter asymmetries, such as the lepton-based ones. The asymmetries are quoted at several 
different stages of the analysis. The easiest number to compare with theory is the "unfolded", or "pro- 
duction level" asymmetry, where the collaborations have processed the measured asymmetry to remove 
background contamination, the effects of analysis cuts and of the detector. Both CDF [6, 7] (9.4 fb _1 ) and 
D0 [8] (5.4 fb _1 ) present their inclusive semileptonic ti asymmetry results at this level. The average of 
the two measurements is 



A ti = 0.174 ±0.038, (1) 

which is significantly larger than the SM prediction, 

= 0.088 ±0.006, (2) 

obtained from next-to-leading order (NLO) QCD and including the leading electroweak (EW) contribu- 
tions [1, 2, 9-1 1]. The SM prediction is derived taking the leading order total cross section in the denomi- 
nator of the asymmetry - a conservative approach - and the error on the SM prediction has been estimated 
by varying renormalization and factorization scales (see e.g.Ref. [12]). 

Both CDF and D0 have also measured the dependence of the asymmetry on the mass and rapidity of 
the tt system. For the m t i dependence, CDF [6, 7, 13] finds, after unfolding: 

A 1 ^ = A t i (m ti < 450 GeV) = 0.084 ± 0.053 , 

A^ gh = %(m tt ->450GeV) = 0.295 ± 0.066 . (3) 

The asymmetry in the high-m t j bin is particularly striking given that the SM prediction, including EW 
corrections [1, 2, 9-1 1], is much lower: 

(4f h) SM s . 129 +o.oo8_ (4) 
For the same quantity, D0 has only reported a measurement without unfolding ("reconstruction level"): 
(A^ gh ) = 0.115 ±0.060. (5) 

V u / reco 



We note that the D0 value for A^ eh is consistent with the CDF value at the reconstruction level, (A^ 
'Flavor- violating, t-channel new physics mechanisms to shift A t i can be Fierz rearranged into a form where this is true 
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0.198 ± 0.043, which suggests that upon unfolding the value obtained by D0 would be larger than the SM 
expectation. Putting it differently, assuming the same unfolding factor between CDF and D0, we find that 
both measurements prefer a rather large value for A^ gh , but with D0 being more consistent with the SM 
prediction. 

The other class of forward-backward asymmetric observables is the lepton-based asymmetries. From 
the same selected events used to measure Ag, both Tevatron experiments have also measured the single, 
Ai, and dilepton, An, asymmetries [13, 8, 14, 15]. The results, given at the unfolded level, 2 are: 

(^) CDF = 0.066 ±0.025, (^h)cdf = 0.42 ± 0.15 ± 0.050 , (6) 

(At) D0 = 0.152 ±0.04, {A U ) D0 = 0.053 ±0.084. (7) 

The SM predictions for the lep tonic asymmetries, as reported by the experimental collaborations, are: 

(A)cdf = °- 016 > (^h)cdf = 0.060 ±0.01, (8) 

(^)d0 = 0.021 ± 0.001 , (A h )d0 = 0.047 ± 0.001 . (9) 

Being a proton-proton collider, the LHC is not sensitive to A t ^. However, the LHC can probe a related 
observable - the charge asymmetry in ti production, Ac- Measurements of Ac at = 7 TeV have 
been reported by both ATLAS [16] (for 4.7 ftr 1 ) and CMS [17, 18] (for 5.0 fb _1 ). The expected SM 
asymmetry [1,2, 11], A C M = 0.0115 ± 0.0006, is much smaller than the Tevatron's asymmetries, due 
to the domination of the gluon-gluon production channel, which is symmetric. The LHC measurements 
so far are consistent with the SM value, but the size of the statistical and systematic uncertainties is such 
that one cannot yet exclude the consistency with the anomalous Tevatron A t i result. It is important to 
emphasize that even within the SM the Tevatron and LHC observables differ in nature. In particular, the 
dominant ti production mechanism and the kinematical reaches available to the top quarks are clearly very 
different at the two colliders; the Tevatron collides charge-asymmetric beams and top quark production 
is dominated by quark-antiquark annihilation, while, at the LHC, collisions are charge symmetric and 
top pair production is driven by gluon-gluon collisions. Furthermore, non-SM dynamics can naturally 
induce a large deviation for the forward-backward asymmetry at the Tevatron without affecting the charge 
asymmetry at the LHC [12, 19, 20]. 

Given that the LHC probes a different observable, we turn our attention back to the Tevatron. The 
discrepancy between the SM predictions and the measured asymmetries at the Tevatron could be due to an 
unknown QCD effect, or an unidentified experimental bias. Alternatively, it might be a hint of dynamics 
beyond the SM (for a review, see e.g.Ref. [21]). Either way, the current situation is not satisfying, and 
the main goal of this paper is to investigate what other information can be used to gain more insight. 
Specifically, we propose a correlation between A t i and daughter asymmetries that are experimentally easy 
to measure and also under theoretical control. This correlation can then be used to distinguish the SM 
from more exotic explanations of A t ^. 

Our basic idea is simple, at least in principle. In the SM, the lepton-based asymmetry in ti events 
is completely determined by the ti asymmetry, meaning for a given A t i one can use top quark decay 
kinematics to predict A\. Radiation originating from the top quark decay products alters the kinematics 



2 Note that the CDF collaboration has not yet published an unfolded result for Ai, the one we state here is given for events at 
the "background-subtracted level". 
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and blurs the relationship between A\ and A t i, however this effect is suppressed by the narrow width of 
the top quark. This relationship is true for the inclusive asymmetries, but also differentially - taking the 
asymmetries with respect to a kinematic variable x; in each bin of x, the lepton asymmetry can be fixed 
knowing A t i in that bin, such that Ai (x) [A t i (x)] traces a calculable curve as x is varied. 

However, once we move beyond the SM, Ai and A t i are generically independent. At high m t i, A[ 
is indeed driven by the top quark kinematics and polarization [22-28], however, near the ti threshold 
A[ is set by the initial-state quark polarization rather than anything related to the top quarks [29]. Thus, 
given some observable x that interpolates between the threshold and high-m t j regions (lepton p?, Ht, 
etc.), the curve Ai(x) versus A t i(x) will be different for models beyond the SM. Our proposal is to use 
x = Pt ( tne lepton pt) and to simultaneously measure A\ (p l T ) and A t i (p l T ), to verify whether the curve 
A[ (p l T ) [A t i (p l T )] is in agreement with the SM. We choose the lepton px as our kinematic variable because 
it is experimentally clean and easy to reconstruct. 

We begin our study of this correlation in Sec. 2, working at the partem level and without cuts to demon- 
strate the basic idea. In Sec. 2.1 we provide several checks that suggest that the correlation is indeed robust 
and therefore more sensitive to new physics contributions to the asymmetries. In detail, we consider three 
types of effects: (1) mis-modeling of the of the tt transverse momenta, (2) scale dependence of the differ- 
ential asymmetries and (3) radiation in the top decay. We verify that, while some of these effects influence 
the overall normalization of the asymmetries, affecting the agreement between theory and experiment, the 
correlation between A\ (p l T ) and A t i (p T ) is unaffected by these deformations. Then, in Sec. 3 we in- 
clude detector acceptance, experimental cuts, parton showering, and top quark reconstruction, to show the 
asymmetry correlation in a realistic hadron-collider environment. Finally, in Sec. 4 we consider several 
simple new physics models and show that the SM correlation is significantly violated in general which 
therefore can potentially lead to much cleaner extraction of a possible non-SM signal. This is followed by 
some discussion regarding the use of reconstruction-free variables (Sec. 5) and our conclusions (Sec. 6). 



2 Idealized case: SM 

To get some intuition regarding our proposal we begin by discussing the idealized SM case where no 
acceptance cuts are included. The differential asymmetry observables are defined as: 

A (J, N AY ti >0 (P 1 T ) ~ ^Ay tt -<o (Pt) nm 
N AY tl >o [Pt) + N AY tt <o [Pt) 

N Yl >o (Pt) + Ny l<0 (p l T ) 

where N&y->o (Ny^o) and N^y tT <o (Ny{<o) are the number of events with AY t ^ {Y{} greater or less than 
zero. 3 We also study the cumulative distributions - the asymmetry for all events with lepton pt above a 
given threshold, obtainable by integrating the numerator and denominator of the differential distributions, 



i Yi is defined as Qi - iji, such that a backwards-moving electron is the same as a forwards-moving positron. 
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then taking the ratio: 



oo 



(N AYti>0 (p l T ) ~ N AYtI <o (Pr)) 



A-tt (PT,cut) 



PT.cut 



(12) 



oo 



(N AYti>0 (p l T ) + N AYti<0 (Pt)) 



PT.cut 



OO 




Al (PT,cut) 



PT.cut 



(13) 



oo 




PT.cut 



The differential distributions contain the physics we want to study - the correlation between A\ and A t j, 
but they are difficult to measure experimentally owing to the limited top quark sample size. Meanwhile, 
cumulative distributions are more tractable experimentally, but the integration over multiple bins dilutes 
the correlation between A t i and A\. We present both types of distributions for the idealized SM case to 
show the similarities and differences. 

We are interested in the lepton asymmetry in the lab frame, as well as the lepton asymmetry after 
boosting to a frame where the ti system has no longitudinal momentum. The lepton kinematics, which 
encode the asymmetry inherited from the top quarks, get smeared under motion of the ti system, hence 
boosting back leads to a larger Ai. The boost only effects the leptonic asymmetry, as A t i is defined in 
terms of a rapidity difference and is manifestly invariant under longitudinal boosts. Unless otherwise 
specified, we will use the generic A\ for the lab frame lepton asymmetry A] ab , and use 
specifically to the lepton asymmetry in the boosted frame. 

The primary tools for our study are the NLO Monte Carlo generators MCfm (v6.3) [30,31], and 
POWHEG (here run in the hardest-emission generator mode) using the heavy quark production routines [32- 
35]. For the idealized SM case, all results were generated using the MSTW2008NLO [36] parton distri- 
bution functions and with factorization and renormalization scales set to /iR = /^f = Q = \f™% + (Pt) 2 - 
Spin correlations between the top (antitop) quark and its corresponding decay products are maintained in 
both codes. 4 

The distributions A t i{p l T ), Ai (p l T ), A t ^(p l Tcnt ) and Ai (p l Tcut ) for the ideal case are shown below 
in Fig. 1. The darker lines show the MCfm results, while the results obtained with POWHEG are shown 
in the lighter shaded lines. The green curves show A\ (p l T ) in the lab frame, while the red curves show 
^boost (j,i T ^ the lepton asymmetry in the Y tf = frame. We find that the two NLO Monte Carlo (MC) 
generators are in reasonable agreement. The qualitative behavior of the curves can be understood as 
follows: beginning with the leptonic asymmetries, near threshold A[ (p l T ) is sensitive to the polarization 
of the incoming quark, which is small due to the vector-like nature of QCD. Hence we expect Ai to be 
near zero [29]. In the other extreme limit, when the lepton's px is very large it has to come from a boosted 



4 Other relevant choices of Monte Carlo generation parameters are m t = 173.0 GeV, T t = 1.31 GeV and nib = 5.0 GeV. 
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Figure 1: The dependence of top quark and leptonic asymmetries on the lepton pt- The left panel shows the 
differential distributions A t f(p l T ), Ai (p T ) for the SM ideal case, while the cumulative asymmetries A t ^ (p l T cut ) , 
A\ {p l T cut ) are shown in the right hand panel. The darker lines depict the MCFM results, while the results obtained 
with POWHEG are shown in the lighter shaded lines. The green curves show Ai in the lab frame, while the red curves 
show Ai after boosting to a frame where the ft system has no longitudinal momentum. This boost only affects the 
leptonic asymmetries. 



top quark, and therefore the lepton-based asymmetry should asymptote to the corresponding value of A t j, 
keeping in mind that within the SM no net polarization is expected for the top quarks in it events. This is 
consistent with the lepton-based asymmetry curves shown in the plot. 

The behavior of A t i (p l T ) can be understood from the lepton pr spectrum in it events and the depen- 
dence of A t i on m t i. The asymmetry A t i (m t i) is a monotonically increasing function of m t j (see Fig. 2 
where for completeness A t i(m t i) and A\ (m t ^) are presented) [5,37,38], however m t j is only weakly 
correlated with p l T (at least up to p l T < 100 GeV). The correlation, at leading order (LO), between these 
variables is shown explicitly in Fig. 3. The weakness of the m t i — p l T correlation implies that the lower 
p l T bins are populated almost equally by a wide range of invariant masses and hence we do not expect a 
significant rise in A t i (p T ) as the lepton p-r is varied. 5 

Higher p l T events are more correlated with large m t j however the large lepton transverse momenta, 
which naively would lead to a large A t i, are forcing those energetic events to be central. As shown in 
Fig. 3 events with p l T > 150 GeV have m t i > 550 GeV, however, with a top quark rapidity difference 
below 0.7 (as shown in Fig. 3 visualizing the correlation between AY^ and p l T ). As central events tend 
to have lower A t i, we actually expect the overall value of A t i and A\ to be below their nominal value 
expected based on the inclusive and high invariant mass values prior to lepton transverse momenta cuts 



5 A similar argument for the flatness of A t i (p't) can be made using the AY t j — p l T correlation (at LO) depicted in the right 
panel of Fig. 3. The differential asymmetry A t j (\AY t i\) rises monotonically/linearly for increasing absolute values of AY t j, 
and from Fig. 3 we see that each p l T bin picks events with a variety of Alfj{ values. The asymmetry in a given p l T bin is thus 
the average of small A t i at low |AYjj| with large A t $ at high |AY t f|. As p l T is increased, the sampling across a range of AYjj 
stays, but the cross section decreases. However, as the change in cross section cancels out in the ratio defining the asymmetry (a 
normalization effect), A t i remains flat. 



7 



(Fig. 2). This is consistent with the distributions shown in Fig. 1. 

We close this subsection by presenting the same result but this time on the Ai — A t i plane where each 
point in the left panel of Fig. 4 corresponds to a different p l T bin, which is nothing but the A\ {p l T )[A t i (p l T )] 
curve mentioned in the introduction. Given the one-dimensional asymmetry spectra explained above, we 
find that the distribution of points is consistent with a nearly vertical line, and both A t i{p l T ) and A\ (p l T ) 
span rather small values. For completeness, we also show, in the right panel of Fig. 4, the asymmetry 
correlation as a function of m t i, even though this is not the main focus of this work. 
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Figure 2: The differential asymmetries A t j (m t f ) , A\ (m t j) for the SM ideal case are shown in the left panel. The 
lepton pt distribution (at LO) is shown in the right panel, normalized to the total cross section. 
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Figure 3: The correlation between the lepton px and the kinematics of the tt system (at LO). The correlation 
between p l T and m t j is shown on the left, and the correlation between p l T and AY t ^ is shown on the right. 
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Figure 4: The correlation between the leptonic and tt asymmetries displayed in the Ai versus A t i plane. No cuts 
have been imposed, and all quantities have been calculated at parton level. In the left panel we show the curve traced 
out as p l T is increased from GeV (bottom point) to 100 GeV (top point) in intervals of 20 GeV - the p l T intervals 
for some of the points are indicated in brackets on the plot. Similarly, in the right panel we present the curve traced 
out as m t i is increased from the threshold value (bottom left) to 1 TeV (top right) in intervals of 50 GeV. 



2.1 Robustness tests 

Having understood the interplay between the SM A t i and A\ in an idealized scenario, we now study how 
various effects impact the Ai (p l T ) — A t i{p l T ) correlations. We continue to work at parton level in this 
section and will not include cuts. Acceptance cuts and more realistic jet description will be studied in later 
sections. 

The first effect we consider is radiation in decay. The b (b) from the leptonic top (or antitop) quark 
can radiate gluons, changing the kinematics and correlations among the top decay products. 6 Radiation 
in decay will obviously not change the top quark asymmetry, but it may impact how A t i is passed on to 
A[, e.g.through the analysis selection cuts. Though radiation in decay occurs at the same order in a s as 
the processes that contribute to A t i, it is further suppressed by the width of the top quark. Recent analytic 
work on the effects of radiation in decay can be found in Ref. [39], and radiation in decay has been 
incorporated into the latest version of MCfm [31]. By comparing the study of Fig. 1 between different 
versions of MCfm, we can study the size of the effect of radiation in decay on our observables A t i (p l T ) 
and A\ (p T ). As can be seen from Fig. 5, the results with and without radiation in the decay are nearly 
identical, indicating that the LO treatment of the top quark decay products (but with spin correlations 
intact) is sufficient to predict the correlation between A t i and A\. As the top quark pr is ambiguous once 
radiation in decay is allowed, we perform this cross-check using a fix-scale choice of Q 2 = m 2 . 

The second test of robustness concerns the pr of the tt system, Pxtt- As * s we ^ known [3], A t i 
strongly depends on p T t j as it controls the level of real emission in the event. Therefore, among other 
effects, mis-modeling of acceptance cuts or biases in the measurement could lead to a change in the 

6 The hadronic W's decay products can also radiate, of course, but since the W is both a color singlet and narrow, this has 
less influence on the lepton kinematics than the radiation from the lepton's sibling b (b) quark. 
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Figure 5: A comparison of MCfm results with and without radiation in the top quark decay. The solid line shows 
the differential distributions A { j (p T ) and Ai (p l T ) (lab frame) as a function of the lepton pt in the most recent 
version (v6.3), while the dashed line shows the same observable calculated with the previous version. The most 
recent version includes NLO effects in decays, while the previous iteration contained all decay spin correlations, but 
no NLO effects in the decays. 



overall normalization of the resulting inclusive and differential value of A t j. Thus, it is important to verify 
whether the A[ — A t i correlation is sensitive to p T t j. Notice that by insisting on large p l T , we are forcing the 
ti+X system into two possible kinematic configurations: (1) the top and antitop quark move in the same 
direction, recoiling against hard initial-state radiation, or (2) the top and antitop quark are back-to-back 
and both are central. In the first case large lepton px are possible since some of the initial and large p T t i is 
inherited by the lepton. In the second configuration, the lepton inherits large pr from the individual top or 
antitop quark, rather than the ti pair. As a result, this configuration is characterized by low p T t i, but large 
m t f. Given the contributions come from such different kinematic regimes, one may worry that our results 
in Fig. 1 come from a delicate cancellation between different effects. Any mis-modeling of or bias in p T t i 
would disrupt such a cancellation and destabilize the correlation. One way to see whether a cancellation is 
occurring is to divide the events into different bins of p T t j and check the A[ (p l T ) — A t i (p l T ) correlation in 
each bin. Again, computing ti production at NLO, we have performed this test and the results are shown 
below in Fig. 6. Quantifying the correlation by the ratio A[ (p l T )/A t ^ (p l T ), foip Tt ^ < 20 GeV, we see it is 
the same as the correlation in events with p T t i > 20 GeV - a strong indication that the A[ (p l T ) — A t i (p l T ) 
correlation does not come from a cancellation of competing effects. Based on our complete NLO analysis, 
we can therefore conclude it is stable against p T t i mis-modeling. 

The final avenue we explore within the SM ideal case is scale variation. Because A t i vanishes at O(a^), 
NLO calculations of differential ti properties only result in a leading order prediction for the asymmetry 
A t i. A more accurate determination of the asymmetry would involve understanding ti, differentially, at 
NNLO. In the absence of this NNLO calculation (for recent progress on this, see e.g.Refs. [40, 41]), one 
estimate of our ignorance regarding higher-order corrections is to vary the scale used in the A t i calculation 
by a factor of two. While we expect the absolute values of A t i and Ai to change as the scale is varied, 
much of the variation is carried in the scale where a s is evaluated, which cancels out if we take the 
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ratio of asymmetries. A stable ratio A\j A t i under scale variation is therefore a sign that the correlation 
pointed out here is robust. In Fig. 7 we show the differential distributions A t i (p l T ), A\ (p l T ) and the ratio 
Ai (p 1 t) I A t t (p l T ) for three different scale choices: Q 2 = Qq, Q 2 = 4 x Qq and Q 2 = Qq/4: where 
Qq = m 2 + (pt) 2 . The ratio is indeed very stable, as can clearly be seen from the bottom panel of Fig. 7. 
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Figure 6: The ratio of differential asymmetries Ai (p l T ) / A t i (p l T ) for different ranges of p T t j. The red and green 
curves show these ratios for events with p T t i < 20 GeV and p T t i > 20 GeV, respectively. The blue curve is the 
ratio across all Px.tt values. Most of the cross section resides in the lower Pt.u region, so the errors on the high-p T f 1 
region are larger because of limited statistics. All calculations were performed in the idealized SM scenario, at NLO 
in top quark pair production. 



3 Realistic case: SM 

We have seen that, in an ideal detector, the lepton and ti asymmetries are correlated and follow a robust, 
predictable curve as a function of the lepton pt- We must now show to what extent this correlation remains 
intact in a true hadron collider environment. We proceed in two steps. First, we continue with a parton- 
level analysis, but impose a set of realistic cuts employed but the collaborations in the actual analysis; 
for concreteness we are going to apply the cuts used by the CDF collaboration, the ones used by the D0 
collaboration are in practice very similar and have negligible impact on our final conclusions. Second, to 
further close the gap with the true experimental conditions, we repeat our study including parton shower 
effects and genuine top quark reconstruction. 



3.1 Parton-level analysis including cuts 

Including possible real emission, the parton-level process for ti production at NLO has 7 final-state par- 
ticles: one lepton, one neutrino, and up to 5 jets, two of which originate from b or b quarks. Inspired by 
CDF [7], we impose the following cuts on these objects. We require: 
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Figure 7: Dependence of the asymmetries on the lepton px for three different scale choices, calculated using 
MCfm. The left, middle and right panel show A t j, A[ (lab frame) and the ratio A\l A t i, respectively. Only the 
differential asymmetries are depicted, however the same trend is present at the cumulative level, cf. Eqs. (10)— (1 3). 
These plots show the ideal SM scenario where no cuts have been applied. 

• exactly one lepton of p l T > 20 GeV and |?/| < 1.1 . 

• $ T > 20 GeV, which we take directly from the px of the neutrino. 

• jets to have Pt, jet > 20 GeV and \r/j e t\ < 2.0, and to be formed with the kx algorithm using a jet 
size of R = 0.4. 7 

• in addition to the above requirements, that all b and b jets are restricted to \r] b \ < 1.0 . This treatment 
of bottom quarks means we are effectively tagging both the b and b quarks, whereas CDF and D0 



7 In the CDF analysis the CDFmidpoint cone algorithm is used, not kr- 
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Figure 8: Left: differential asymmetries A t i(p l T ) and A[ (p T ), where all CDF analysis cuts have been imposed. 
The darker shaded lines indicate results obtained via MCFM, while the lighter ones show the POWHEG results. Both 
A t t and Ai were calculated in the lab frame, with Ai depicted in green. The blue lines show A t j calculated using 
reconstructed (anti)top quarks, following the scheme in Ref. [31], rather than using the truth information. The grey 
line shows A-g including cuts, but calculated using the true top quark four-vectors rather than using reconstructed 
objects. Right: the A\ (p l T )[A t i (p T )] post-cuts curve with (red line) and without (blue line) reconstruction compared 
to the A;[A t j] curve obtained in the ideal case (grey). Because of the lepton cut imposed in the CDF analysis, the 
new Ai (p T ) \A t i (p T )] curve starts at 20 GcV rather than zero. The brackets indicate the p l T range corresponding 
to particular points along the Ai (jx T )[A t t (p?)] curve. The grey arrow indicates how the Ai, (jr T ) [A t j [p T )\ curve 
shifts from cut effects, while the blue arrow shows the effect of reconstruction. Note that only A t i is affected by the 
reconstruction, while both asymmetries are decreased by the analysis cuts. 



only require one tagged jet. This small difference in acceptance may lead to a difference in the 
absolute values of our calculated asymmetries compared to experiment, but we do not expect it to 
change our main result. 

• isolation criteria to be satisfied: Ai?j cti j c t > 0.4 and Ai^j ct > 0.4 . 

Keeping the jet and missing energy cuts fixed to the above, we calculate A t i (p T ) and Ai (p T ), just as 
we did in the ideal case. As a first step to evaluate this asymmetry, and in particular to evaluate A t i (p l T ), 
we take the top quark truth information directly from the simulation (no reconstruction). The result is 
shown below in Fig. 8. While the absolute value of the asymmetries is diminished once cuts are imposed 
(see Fig. 1 for a comparison), the trend in Ai, A t i versus p l T is preserved. 

Another step towards a more realistic analysis is to redo our top quark analysis in terms of reconstructed 
objects rather than using the truth partonic information (before the decay). While a completely realistic 
reconstruction will be presented in the next section and will include parton showering effects, we can 
study some reconstruction effects even at the parton level. Specifically, even at the parton level, jets 
are sometimes lost due to acceptance or gained from ISR (initial-state radiation), leading to an incorrect 
reconstruction and a warped A t p To study this effect, we compare the A t i calculated using perfect top 
quark reconstruction with the Ag calculated using a reconstruction where all 5 jets, not just those coming 
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from the top quark decays, are considered. This comparison is shown in Fig. 8 as well. The impact it has 
on A t i is sizable, amounting to a 0(30%) reduction in A t j. We find that both effects, namely including the 
acceptance cuts as well as the top quark reconstruction reduce the resulting asymmetries. This is expected 
since they force the events to be more central, and wrong partonic assignment in realistic reconstruction 
dilutes the asymmetry. 

Finally, we have repeated the scale variation check with the post-cut parton-level events and find the 
same trend as in the ideal cuts case (cf. Fig. 7): while the individual asymmetries A t i (p l T ) and A\ (p l T ) 
shift by 0(25-30%) as the factorization/renormalization scale is varied by a factor of two, the ratio 
Ai (Pt)/Au (Pt) remains stable. 

3.2 Showering and reconstruction effects: POWHEG + PYTHIA and SHERPA's 

CSSHOWER 

Having included cuts in our analysis and observed that the correlation between A t i and A[ is maintained, 
the next step towards reality is to include a parton shower. The radiation, which quarks and gluons emit as 
they lose energy in their evolution from the scale of the hard process to the hadronization scale can show 
up as additional jets in the detector. This spray of energy, and the combinatorial problem it creates in any 
analysis relying on reconstruction, tends to dilute parton-level effects. It is therefore important for us to 
show that our correlation remains visible in the actual environment where the experiments reside. 

One way to study the effects of the parton shower on our observable is to pass NLO POWHEG events 
through PYTHIA [42]. Manipulating the particle-level PYTHIA output into jets via FastJet [43, 44], 9 we 
can then apply the same analysis cuts as in the previous section. The only subtlety is how we handle 
6-tagging. In the simple reconstruction used in the last section, we assumed the identity of both b quarks 
was known, leaving no ambiguity about whether/how to combine the tagged jets with the lepton. In 
this section, we drop this assumption, choosing to select and reconstruct events exactly as in the CDF 
analysis [13] . Specifically, provided an event contains a lepton and missing energy passing criteria shown 
in Sec. 3.1, the event is kept if it contains four or more jets and at least one 6-tag. The leading four jets, 
one of which must be a b jet, are subsequently divided into a hadronic top quark system and one jet to be 
paired with the lepton + $ T system. Each combinatoric possibility in the leptonic W reconstruction and 
the division of jets is tested, and the combination that best reconstructs the top quark masses (the hadronic 
and leptonic) is retained. The differential asymmetries for the reconstructed top quarks (and leptons) are 
shown in Fig. 9 below, and in the A\ (p l T ) [A t i (p l T )] plane on the left panel of Fig. 10. For comparison, we 
also show A t i (p l T ) and A\ (p l T ) at parton level. 

The experiments typically compare asymmetries at the background-subtracted and unfolded levels, the 
analog of comparing our fully showered and "ideal scenario" results. However, as we are dealing with 
Monte Carlo data rather than experimental data, we are able to break down these shifts in multiple stages. 
Being more quantitative (cf. Fig. 8), we have seen that imposing cuts reduces the leptonic asymmetry by 
roughly 50% and A t i by 20% when compared to the "ideal case" values. Reconstruction does not affect 
Ai, but it further reduces A t i by an additional 30%. Working at the parton level, these reductions are 
almost independent of the lepton px, our proxy for the energy of the ti system. Showering and more 

8 Specifically, we use the "improved" reconstruction introduced in Ref. [31]: all non-6 jets are considered in the hadronic W 
and t reconstruction, and the combination (2 or 3 partons) that has invariant mass closest to mw or mt is selected. 

9 We use the kr algorithm with jet size R = 0.4. Jets are identified as b jets if there is a 6 (6) parton within the jet area. 
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Ai (p l T ) at parton level. 



realistic reconstruction have a p T dependent effect; A[ is shifted by 0(1-10%), and A t i is shifted by 
0(20-30%), with the shifts increasing with p l T . 

As a second study of the effects of a parton shower on the A\ — A t i relationship, we analyze the 
asymmetry correlation using Sherpa (v 1.4.0) [45]. Despite the fact that Sherpa is a LO matrix element 
generator, the addition of a parton shower technique that appropriately includes color-coherence effects 
will generate a forward-backward top quark asymmetry (see Ref. [46]), so it is interesting to see what 
degree it retains the A[ — A t i correlation. To generate events, the LO matrix elements, including the 
decays of the top quarks, are showered according to Sherpa's color-coherent showering (CSSHOWER) 
description [47,48]. In the infrared limit of QCD, this description correctly accounts for ISR effects, 
intermediate top quark radiation and multiple emissions from the final-state b quarks and decay partons. 10 
While qualitative agreement in A t i between the LO matrix element plus parton shower (LO + PS) results 
and the complete NLO results has been demonstrated [46] , exploring Ai — A t i tests the adequacy of the 
LO + PS calculation on a detailed and differential level. 

The Sherpa events are reconstructed in similar, but slightly different manner than we used in the 
POWHEG + PYTHIA case. In particular, the charge of the b jets is assumed to be known, leading to perfect 
assignment of the b (b) quarks to the top (antitop) quark. 1 1 After combining the lepton and missing energy 
with the correct b jet, the top quark objects are reconstructed by testing all possible partitions of jets - 
even beyond the 4 jets required for event selection - and keeping the combination which yields the smallest 

l0 For the hard process generation, we use CTEQ6L1 PDFs [49] and an mx,t-hke scale choice (m\ t = mf + (pt) 2 ) 
resulting from utilizing the default scale setting prescription applied in SHERPA. The top quark decays are incorporated at full 
matrix element level, i.e. pre serving spin correlations and full width effects beyond the narrow-width approximation, with the 
only requirement of producing intermediate top quark states (m t = 173.2 GeV). 

11 We obtain this information from a kr jet algorithm where we systematically keep track of the &-parton flavours. 
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Figure 10: Left: the p l T dependence of the leptonic and tt asymmetries for POWHEG events that have been showered 
and hadronized with PYTHIA. The post-shower Ai (p l T ) [A t t (p l T )] curve is shown in dark green and is overlayed on 
top of the result with cuts and reconstruction done at parton level (red). The blue line shows the post-shower 
result in events where the total Pxtt is restricted to Pt tt < 10 GeV. Right: the analogous curves for events 
generated with Sherpa's CSshower. The red curve shows Ai (p l T )[A t i (p 1 t)) witri tne CSshower truncated at 
the first emission, the blue curve shows the fully showered result for events restricted to p T f j < 10 GeV, and 
the green curves depict the inclusive, fully showered SHERPA results. The difference between the green curves is 
the reconstruction algorithm. In the dark green curve, as described in the text, the number of jets included in the 
reconstruction algorithm is allowed to float, while in the light green curve the number is restricted to four, giving 
a straightforward, unambiguous reconstruction out of the hardest b and b jet and the two leading light-flavour jets. 
The low-p T t i result uses the latter reconstruction algorithm. 



summed-up mass deviation, which we take as | m^ eudo _ t -m t \ + \ m^ udo _ t - m t \ + \/2 | m do . t Jets -m w \. 
By including more jets, the combinatorial issues are bigger for this reconstruction method, however it will 
capture situations the CDF method cannot, i.e. where the hadronic top quark radiates and manifests itself 
actually a 4-jet system, or the jet system of the leptonic top quark possesses an extra jet beyond the b jet. 
The leptonic and ti asymmetries derived from the reconstructed SHERPA events are shown in the right panel 
of Fig. 10. We also show the asymmetries obtained from running Sherpa's CSshower in "one-emission" 
mode, where showering is terminated after no more than one emission has occurred. The one-emission 
mode allows us to see the impact of the multi-parton emissions and the (enhanced) combinatorial headache 
coming with them. 

Inspecting Fig. 10, while the shapes of the differential asymmetry curves at parton (or "one-emission") 
level agree reasonably, the results after the full shower are quite different. The slope of the POWHEG + 
PYTHIA curve acquires a slight tilt, as the later/secondary shower emissions drive the slope of A t i{p l T ) 
and Ai (p l T ) slightly in opposite directions (see Fig. 9). For the SHERPA result, higher/secondary emissions 
drive A t i and Ai in the same directions as POWHEG + PYTHIA (the A\ is slightly drecreased, while the A t i 
is enhanced, 12 especially at high p l T ), however the change in slope is much more dramatic: A[ is decreased 
nearly to zero, or even negative, until p l T >!00 GeV. Massaging the SHERPA top reconstruction algorithm 



A larger enhancement from Sherpa's CSSHOWER is expected [46] because subsequent emissions in the shower are still 
carried out in a color-coherent manner. 
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does not fix the discrepancy - the simplified reconstruction used in the lower curve of the Sherpa panel 
in Fig. 10 is closer to the method used for POWHEG, yet the A\ (p l T )[A^ (p l T )] curve is just as disparate. 

The \ow-p T t i curves in Fig. 10 offer some insight into why the post-shower results disagree. The trans- 
verse momentum of the ti system separates the so-called "Sudakov region" - where A t i is positive - from 
the "hard-p-r region" where the entire ti system recoils against additional radiation and the asymmetry is 
negative. By imposing a low cut on p T t i we can focus on the effects that drive the positive asymmetry and 
are generated by the virtual corrections and multiple soft gluon emission. The majority of the cross section 
resides in this low-p^ region, so it is important to understand and study the asymmetry here. As we can 
see, the low Pxtt curves closely follow the ideal NLO case (cf. the discussion in Sec. 2), despite the fact 
that all virtual pieces beyond the soft and collinear approximation, i.e.beyond the lowest order Sudakov 
description are absent in the Sherpa CS SHOWER description. In other words, for regions of phase space 
close to Born kinematics, the LO + PS and NLO calculations are in relative agreement. Note that the 
denominator in the Sherpa asymmetry calculation is the lowest order cross section and is smaller than the 
NLO cross section used for POWHEG. 

However, turning to the inclusive case, the calculation becomes sensitive to the whole p T t i region and 
how it is modeled. As discussed in [46], in color-coherent parton showers like the CSSHOWER, the real- 
emission effect is overestimated because leading- Nc color factors, rather than the smaller, NLO-correct 
factors, are utilized and the use of the eikonal limit/dipole radiation functions is extended beyond its/their 
reliable range. The enhanced emission leads to a stronger reduction of the already smaller virtual effects as 
soon as we go away from the low p T t i limit. The shower Sudakov effect (the enhanced asymmetry in the 
low-p T 1 1 region) is weaker than the enhancement in the full NLO calculation, so the negative contributions 
from emission win out and drag the asymmetries down. So, although the Sherpa CSSHOWER qualitatively 
generates an A t i{p Tt i) similarly to POWHEG + Pythia, it cannot reproduce it in the details, and the 
details are important to get the Ai (p l T ) [A t i (p l T )] correlation right. One reason the lepton asymmetry may 
be more sensitive to the shower details may be the tight rapidity cut, \rj\ < 1.1, whereas |?7j c t| < 2.0. 
Being so central, the leptons can easily be nudged from forwards-moving to backwards-moving (or vice 
versa) or driven out of (into) the acceptance range, making A[ more susceptible to later shower emissions. 
We emphasize that the instability shown in Fig. 10 is not SHERPA-specific, any LO + PS calculator will 
struggle to capture the details of the A\ [p l T )[A t i (p l T )] correlation for the reasons mentioned above. 

To summarize, we find that the qualitative features of the SM A t i (p l T ) — A\ (p l T ) correlation are main- 
tained throughout all levels of event and analysis complexity. The absolute values of the asymmetries 
and their relation to each other do shift depending on the stage of the analysis, but the slope of the 
A[ (p l T )[A t i (p l T )] curve is maintained. Importantly, to see the robustness of this result, one must use 
calculational tools that are exact at NLO (such as POWHEG + PYTHIA). Lowest order matrix element cal- 
culators supplemented with a parton shower respecting color coherence may be sufficient to describe gross 
features of A t i and do offer valuable insight to the physics behind the asymmetry, but these tools lack some 
of the physics necessary to capture detailed effects like the correlation between differential asymmetries. 
While the sensitivity of differential correlations shown here should serve as a warning label on calcula- 
tions done with LO + PS accuracy, as we will show in detail in the next section, the Ai (p l T )[A t ^ (p l T )] 
curves in various benchmark BSM scenarios are dramatically different from the SM curve. Therefore, 
small shifts (in absolute value) in the asymmetries originating from cut, reconstruction or shower effects 
will not seriously hinder the discriminating power of the A\ (p l T ) [A t i (p l T )] correlation. 
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4 Asymmetries beyond the standard model 

As we have discussed, in the SM the shape of A\ (p T )[A t ^ (p T )] can be simply understood as a conse- 
quence of the fact that tt production in QCD is unpolarized. In particular, the same number of left- and 
right-handed top quarks is produced, and that equal amounts of ti pairs are produced in collisions of left- 
and right-handed quarks. The latter fact ensures that the lepton asymmetry vanishes at the tt production 
threshold [29]. On the other hand, the fact that the top quarks have no overall polarization ensures that, 
for top quarks produced with a significant momentum, the lepton direction is determined by the kinemat- 
ics of the boosted top quark decay. These simple arguments explain the behavior of A\ (p l T )/ A t i (p l T ) in 
the SM. At small p l T the lepton asymmetry is dominated by the near-threshold sample, and one expects 
M (Pt) ^ A t t(p l T ). As p l T grows the lepton direction becomes more and more correlated with that of 
the parent top quark, therefore one expects A\ (p l T ) to approach the A t i(p l T ) curve from below. This is 
indeed what comes out of the MC simulations we presented in Sees. 2 and 3. Moreover, we found that 
the top quark asymmetry is approximately independent of p l T due to an interplay between two effects: the 
increase of A t i as a function of m t i (correlated with p l T ), and the increase of A t i as a function of AY t i 
(anti-correlated with p T ). 

These expectations can be grossly violated in models beyond the SM. New physics may introduce two 
effects that could potentially distort the shape of A\ (p l T ) [A^ (p l T )] away from the SM prediction. One 
is that BSM models addressing the anomalous tt asymmetry at the Tevatron may lead to a very different 
m t i dependence of the asymmetry. Due to the correlation between p l T and m t i, in the presence of new 
physics there is no reason for A t i{p l T ) to remain approximately constant as in the SM. In fact, as we 
will see below, A t i (p l T ) can either increase or decrease, and the strength of the effect is strongly model 
dependent. The other important effect is polarization. Models addressing the anomalous top quark asym- 
metry measured at the Tevatron introduce new particles with different couplings to left- and right-handed 
quarks. This typically leads to polarization both in production and in the final state. These polarization 
effects may be observable, as pointed out many times in the literature [50, 28, 24, 51-53, 27]. In particular, 
the inclusive [25] and threshold [29] lepton asymmetries are important observables to test the SM and 
discriminate between different models of new physics predicting the same top quark asymmetry. Here we 
point out that studying the p l T dependence of A t i (p T ) and Ai (p l T ) provides further discriminating power. 

To illustrate it, we study A t ^(p l T ) and A\ (p l T ) in BSM models that lead to an enhanced top quark 
asymmetry in agreement with the Tevatron observations. We focus on the so-called s-channel axigluon. 
In this model one introduces a color-octet field G a with mass mg and the couplings to quarks that are 
assumed to be flavor diagonal but otherwise arbitrary: 

C C g L)i Ht Gl T a qi + 9R:l qi <f G«T a ($, (14) 

where and qf denote left-handed doublet and right-handed singlet quarks, respectively. If the axigluon 
couplings to the top and light quarks are chiral, ql 9r, then interference of the s-channel axigluon 
exchange with the tree-level QCD gluon exchange results in a non-zero forward-backward asymmetry 
of the qq — > ti production process. There exist regions in the parameter space of the axigluon mass and 
couplings leading to an 0(10%) additional contribution to the inclusive top quark asymmetry, bringing the 
theoretical prediction to a good agreement with the central value observed by CDF and D0, while keeping 
the total ti cross section within the experimental bounds. To avoid the bounds from dijet production and 
ti resonance searches, the axigluon should have a significant width, Fq > 0.2 tug- 
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In the following we study in detail the predictions concerning A t i (p l T ) and A\ (p l T ) for a representative 
set of axigluon benchmark models. One class has a relatively light and wide axigluon, mc = 200 GeV, 
Tq = 50 GeV GeV, with flavor universal couplings to the SM quarks. For this mass, we pick three 
benchmarks with left-handed, right-handed, and axial axigluon couplings to the quarks, each predicting 
the same AA ft - = 0.12 contribution to the inclusive top quark asymmetry (on top of the SM contribution 
of 0(9%)), but each predicting a different lepton asymmetry. The couplings, the axigluon width, and the 
computed lepton asymmetries take the following values: 

(L) g Rti = 0, g Lti = 0.8 g s : AA l = -0.07, 
(R) gR>i = 0.8 g 8 , g w . = : AAi = 0.18, 

(A) g R>i = 0Ag s , g Ul = -0Ag s : AAi = 0.05, (15) 

where g s is the QCD coupling. For those benchmarks the axigluon decay width into qq pairs is of order a 
GeV, thus the larger width we assumed must be explained by exotic decay channels (see e.g. [54]). 

The other class of benchmarks has a relatively heavy axigluon, = 1.5 Te V, and flavor non-universal 
couplings (see [55-63] for some theoretical construction of heavy axigluon models). Again we pick 3 
benchmarks with left-handed, right-handed, and axial axigluon couplings: 

(L) g Liq = -1.3 g s , g R , q = 0, g L .t = 6g s , 9R,t = : AA l = -0.01, Tq = 970 GeV , 

(R) g L , q = 0, g R , q = -l.lg„ g L ,t = 0, g R ,t = 6g s : AA t = 0.14, T G = 460 GeV (16) 

(A) g L>g = 0.6 g s , g R ^ q = -0.6 g L>t = -3g s , g R>t = 3 g s : AA t = 0.06, T G = 350 GeV 

where above q stands for doublet and up-type singlet quarks of the first and the second generation (the 
couplings to down-type singlet quarks are assumed to vanish). In this case the axigluon width is assumed 
to be dominated by 2-body decays into the SM quarks, as follows from the coupling values listed above 
Much as for the light axigluon, each of the benchmarks in Eq. (16) predicts AA t i = 0.12. 

In Fig. 11, we plot the p l T dependence of the top and lepton asymmetries. To compute the asymme- 
tries, we simulated semilep tonic ally decaying tt events using Madgraph 5 [64] with a custom user-defined 
model describing the extension of the SM with the axigluon coupled to quarks as in Eq. (14). The asym- 
metries were computed using the parton level input without taking into account showering, hadronization, 
detector, or reconstruction effects. The samples were divided according to p l T into 6 bins with the lower 
bin limits at 0, 20, 40, 60, 100, 150 GeV, and the points in Fig. 11 refer to the asymmetries in these bins. 
As our simulations are tree-level and do not include the 1-loop SM contribution, we simply add to the 
results in each bin the SM asymmetries in that bin estimated by POWHEG. 13 

As can be seen in Fig. 1 1 , the differences between the bechmarks become more pronounced when 
p l T dependence is explored. The striking observation is that the slope of the curve traced in the A^- 
Ai plane can be completely different than in the SM. For the light axigluon case, the (R) benchmark 
predicts an approximately constant A% (p l T ) and increasing A t i (p l T ), while for the (L) benchmark A[ (p l T ) 
is increasing and A t i (p l T ) is decreasing. In both of these cases the shape of the curve can be qualitatively 
understood. At low p l T the lepton asymmetry is close to its thresholds value and approaches A t i (p l T ) at 



"Both the SM and BSM asymmetries are dominated by interference of the tree-level QCD amplitude: with the 1-loop ampli- 
tude in the former case, and with the tree-level s-channel axigluon exchange in the latter, so it is reasonable to assume that the 2 
effects add up. 
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Figure 11: Left: the At (p l T )[A t i (p l T )) curves for increasing p l T bins in the axigluon benchmarks with m,Q — 
200 GeV with left-handed (blue), right-handed (green), and axial (red) couplings, as specified in Eq. (15). Also 
shown the analogous SM curve (black) for comparison, and a "rescaled SM" curve (grey), where A t i has been 
artificially inflated to 0.18 - the current experimental value. This rescaled curve emphasizes the difference in slope 
between the SM and our benchmark BSM scenarios. Right: the same curves presented for the axigluon benchmarks 
with me = 1-5 TeV as specified in Eq. (16). 



large p l T . The shape of A t i (p l T ) is influenced by top polarization: in the (L/R) case boosted top quarks 
have dominantly left/right helicity and tend to emit the lepton in the opposite/same direction as the top 
quark. For low and moderate p l T this leads to a certain degree of anti-correlation/correlation with m t j, and 
the derivative of A t i (p l T ) reflects the increasing tt asymmetry as function of On the other hand the 
(A) benchmark, where polarization effects are small and completely vanishing at the ti threshold, predicts 
the shape A t i (p l T ) and A\ (p l T ) similar to the SM, except that both lepton and top asymmetries in each bin 
are shifted to larger values. 

Very similar arguments hold for the heavy axigluon benchmarks. One important difference between 
light and heavy axigluon benchmarks is that the latter predict a much steeper dependence of the top 
asymmetry on m t ^. This is the reason why p l T dependence is more pronounced for the heavy axigluon 
benchmarks. In particular, A t i (p l T ) steeply grows at high p l T as a result of the correlation between p l T and 
m t p For the (L) benchmark this leads to a "turnaround" of the curve when anti-correlation between p l T 
and m t i at low p l T due to top polarization turns into correlation at high p l T . 

In summary, p l T dependence of the lepton and top asymmetries offers a handle to discriminate between 
the SM and new physics, and also between different models of new physics predicting the same top asym- 
metry. This test of the SM is to a large extent independent of the overall normalization of the asymmetries. 
If, hypothetically, some yet un-calculated higher-order QCD corrections to the asymmetries happen to be 
much larger than expected, they may shift the curve in the Ai (p l T ) — A t i(p l T ) to higher values without 
changing its shape. On the other hand, the shape of the curve in the A\ - A t i plane is very sensitive to 
polarization effects, which typically are present in BSM models addressing the anomalous top asymmetry 
at the Tevatron. 



20 



5 Discussion: additional reconstruction-independent quantities 



Throughout this paper we have used p T to study the relation between Ai and A t i in various kinematic 
regimes. The lepton pr is a useful variable to consider as a probe since it is simple, clean, and - just like 
Ai - it does not require reconstruction of any complicated objects. However, it is not the only option. 
Along with the properties of the lepton, there are several global observables, such as the total invariant 
mass taken over all visible final-state objects, m V i s , that are also reconstruction-independent. A few other 
examples are Ht, rriT and mT,vis- Here, Ht is defined in the usual way as the scalar sum Ht of all 
identified-object including $ T ; similarly, combining all leptons and jets in the event (even beyond 
those required for selection) into a single "visible" four-vector, p v ; s , the other transverse quantities are 
given as (cf. Ref. [65]): 

m T = E WtA + \$t\\ - (pt^s + Bt) 2 , 

\i=J,jet / 

"4,vis = ™vis + 2 (E T , v is \$t\~ PT.vis • $ T ) (17) 

where 

m lm = Pvis> PT.vis = PT>i> E T, vis = ™vis + PT, vis ( 18 ) 

i=l, jet 

denote the mass, transverse momentum, and transverse energy of the visible system. 

Studies using reconstruction-independent asymmetries and/or variables such as these are useful for a 
couple of reasons. First, reconstruction-independent quantities have no (intrinsic) combinatorial issues, so 
the sensitivity to how these objects are chosen is significantly reduced. Second, as they are inclusive, these 
quantities are less sensitive to the modeling of the various distributions involved in the actual measure- 
ments. As a result, systematic uncertainties are under better control. Third, we can introduce additional 
analysis levels with varying degree of reconstruction-dependence (e.g.reconstruction-dependent asymme- 
try, A t i, versus reconstruction-independent observable, mr), allowing for a number of sanity and closure 
tests on the existing data. For example, if the Tevatron collaborations were to measure the dependence of 
the asymmetries on one of the variables introduced above and find agreement with our studies, that would 
strongly hint at an error in the reconstruction. With no more data-taking at the Tevatron it is essential to 
try out each different angle in forming differential asymmetry measurements as it might be the only way 
to shed more light on the present data versus theory puzzle. 

To illustrate this in more detail, we present in Figure 12 the differential top quark (left panels) and 
lepton (right panels) asymmetries with respect to the final-state variables discussed above. These events 
were generated using POWHEG + Pythia (top panels), and simultaneously with Sherpa (bottom panels) 
following the same procedure as in Section 3.2. For comparison, we also show the differential asymmetries 
with respect to m t f. The qualitative similarity between results in the different generators is evidence for 
the decreased modeling sensitivity mentioned above. 

Focusing on A t i, the asymmetry with respect to Ht does not rise nearly as fast as A t ^(m t A. This 
occurs because Ht is sensitive to extra radiation in the event, while m t i is not. Events that have large 
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Figure 12: Differential forward-backward (left) and lepton-based asymmetries (right) in tt events as a function 
of various transverse observables (and for comparison, as a function of m t j) using CDF-like event selection (see 
Sec. 3.1). The results in the top two panels were obtained by running POWHEG + Pythia, while the bottom panels 
show results from Sherpa (all generated as in Sec. 3.2). 



Ht due to lots of extra radiation will carry a large, negative asymmetry, partially counteracting the large 
positive asymmetry generated by large m t i, large Ht events. Similarly, the design of the selection cuts 
is such that in most cases the visible transverse momentum is balanced out by the missing transverse 
energy vector. This yields an asymmetry dependence on tut that is very close to that of Ht- Finally, we 
use mr,vis rather than m v - ls alone, because mr, v i s contains some missing energy information and thereby 
gives a better approximation to the total invariant mass of the tt+X system. The effect only kicks in 
for larger transverse masses where A t i (mr lV is) starts to follow the m t i dependence of A t j. In the low 
mT,vis region the transverse character of the observable still dominates, leading to an ify-like asymmetry 
dependence. 
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Turning to the dependence of Ai, we find that the differential leptonic asymmetry actually is quite 
insensitive to the choice of the specific transverse variables and is very small in magnitude (\A\ (x) | < 2%), 
potentially too small to be observed at the Tevatron (even without requiring top quark reconstruction). 
However, this argument can be turned around in the sense that the standard model gives predictions that 
by and large are "consistent with zero". Hence any measurement deviating from zero in these observables 
can be considered as a fairly clear signal for new physics. 

It is interesting to explore these "reconstruction-free" distributions both from an experimental point of 
view as well as from the theoretical perspective, where one should study the robustness of this new set 
of observables more carefully. Such a semi-inclusive approach can, in principle, lead to a cleaner set of 
precision observables and, in the future, also turn out to be useful for the LHC experiments. 



6 Conclusions 



The ti asymmetry measured by the CDF and D0 collaborations remains higher than expected in the 
standard model (SM), both when measured inclusively and when binned in kinematic variables such as 
m t i or Y t p In order to determine the origin of the discrepancy, in this paper we have proposed that the 
leptonic asymmetry Ai is a useful probe. Within the SM, A\ is inherited from A t i, with the relationship 
between the two set by the kinematics of the ti system. Thus, by varying the lepton pt - a simple proxy 
for the mother-top quark energy, directionality and spin - one obtains a simple correlation between the 
top quark and lepton-based asymmetries which can be conveniently described by a curve Ai[A t i (p l T )] in 
the Ai — A t i plane. We have verified that this correlation between asymmetries is qualitatively maintained 
through all levels of the analysis, from partem level through the inclusion of showering and reconstruction. 
The correlation is also stable under variations of theory inputs and under potential mis-modeling. 

By studying the full A\ — A t i correlation in data (and their individual p l T dependence), rather than 
just the inclusive asymmetries, we gain discriminating power. The slope of the curve is insensitive to 
the size of the overall normalization of the asymmetry and is found to be rather stable against various 
deformation of inputs and analysis parameters. The slope is, however, not satisfactorily modeled by LO + 
PS tools and requires complete NLO (+ PS) treatment. We have demonstrated the discriminating power 
of the A\ {p l T )[A t i (p l T )] correlation by showing several beyond the SM benchmarks with the same net 
asymmetry but with correlations in the A\ — A t i plane that are all distinct and distinct from the SM 
case. This argument for studying the correlation among asymmetries is not unique to the semileptonic ti 
channel, nor is it unique to the Tevatron. It is important to keep in mind that, even if the central values of 
observables such as Ac at the LHC turn out to match the SM values within uncertainties, the correlation 
of Ac with other asymmetries in the ti system are a worthwhile test of the SM, and can be proven to be 
more stable similar to the case discussed above. We have also briefly considered how the asymmetry can 
be measured via more inclusive, "reconstruction-free" variables. These can be further used in principle to 
control the measurement of top quark precision observables and might provide us with extra sensitivity to 
the presence of non-standard model contributions to the various asymmetries. 
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